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Abstract — Consider a wireless network of transmitter-receiver 
pairs. The transmitters adjust their powers to maintain a 
particular SINR target in the presence of interference from 
neighboring transmitters. In this paper we analyze the optimal 
power vector that may achieve this target in the presence of 
randomness in the network. Specifically, we start from a regular 
grid of transmitter-receiver pairs and randomly turn-off a finite 
fraction of them. We apply concepts from random matrix theory 
to evaluate the asymptotic mean optimal power per link, as 
well as its variance. Our analytical results show remarkable 
agreement with numerically generated networks, not only in one- 
dimensional network arrays but also in two dimensional network 
geometries. Remarkably, we observe that the optimal power in 
random networks does not go to infinity in a continuous fashion 
as in regular grids. Rather, beyond a certain point, no finite 
power solution exists. 

I. Introduction 

Transmitted power is an important resource in wireless 
networks and therefore power control has been crucial since 
the development of legacy networks. For example, the intro- 
duction of efficient power control algorithms (both closed loop 
and open loop), was one of the main improvements third 
generation CDMA-based cellular networks brought about. 
Several algorithms have appeared that provably allow users 
to obtain e.g. a minimum SINR requirement SINRk > jk 
for link k while minimizing the total power or the power per 
user, subject to the feasibility of this solution. |[T|, IS) 

Of course, power optimization remains an important prob- 
lem in emerging and future networks. Ad hoc networks are one 
such class, where substantial effort has been made to analyze 
their behavior, such as connectivity and transport capacity |3|. 
However, the majority of such work has been focused on fixed 
power transmitted and has not analyzed the potential benefits 
of power optimization. Power control together with scheduling 
becomes more important when multiple hops on the network 
are necessary for information to reach its destination. In a 
different setting, the density of WiFi networks is already high 
enough to create significant interference. This will also be- 
come an issue when femto-cells, a recently proposed network 
paradigm, will become massively deployed: due to their close 
proximity, neighboring femto-cells may create interference to 
one another Q . As a result, all above situations are expected to 
benefit tremendously from power optimization. Nevertheless, 
little progress has been made in finding analytic estimates 
of the performance of random interference-Umited networks 
when power control is apphed 0. 
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In this paper we apply tools from random matrix theory, 
(6], fT\ provide an analytic estimate of the optimal power 
performance for a large network in the presence of both 
interference and randomness. We will start from an ordered 
network structure, in particular an equally spaced line (or 
square) of N transmitters, each with a receiver located in its 
neighborhood at a fixed distance. We will then add randomness 
in the network by removing each base-receiver pair with 
probability c. As a result, a network with roughly A^(l — e) 
transceiver pairs will remain, randomly located in the original 
grid. This sparse network is a model for a realistic cellular 
network, which usually tries to mimic a (usually hexagonal) 
symmetric network. It is also a good model for a wireless 
network with intermittent activity, such that at any time only 
a fraction p of the total bases are active. Although relatively 
straightforward, this paradigm carries all the traits of a random 
network and, as we shall see, behave qualitatively different 
from ordered networks. Also, surprisingly these results are 
valid for both one-dimensional and two dimensional networks 
(arrays and square-grids). 

The paper is organized as follows. In the next section we 
define the network model in the absence of any randomness. 
In Section |III] we describe the network with erasures of base- 
receiver pairs and calculate the average minimum power per 
user as well as its variance using results from random matrix 
theory. We also discuss the relevance of these results in both 
one-dimensional and two dimensional geometries. Finally, in 
Section HV] we draw some conclusions, while at the Appendix 
we discuss some details of the proofs. 

II. Model Description 

Assume a network of N base-stations on a lattice. The 
lattice may be one-dimensional with equal spacing £ between 
them, or two-dimensional (see Fig. [U. In the latter case, 
which is obviously the more realistic one, it may have any 
symmetry, e.g. hexagonal or triangular, but for concreteness, 
we will assume that the lattice is square. Each base-station 
is connected via a link to a single user located at a distance 
5 from its corresponding base-station. The (average) channel 
coefficient between base-station i and user j is given by 

where nij is the lattice vector of integers corresponding to 
base-station i. For simplicity, we have normalized the channel 
gains to unity for i = j. We also define s — S/£. For d- 
dimensional lattices (d — 1, 2), it can be expressed in terms of 



integer linear combinations of the d basis-vectors of the lattice 



m, 



d 

E 

k=l 



(2) 



for integer < Wfc < L, where N — Since the integers 
ruk fully specify the position mi, we will drop the index i 
when possible. For simplicity we assume periodic boundary 
conditions on the lattice, hence mk = mk + L. This means that 
the distance between any two points is taken as their minimum 
distance on a toroidal geometry (or a circular geometry for 
one dimension). This has no effect on neighboring sites or for 
large sizes, and, in addition, it enforces the circulant property 
of the matrix gij . Also, a > 2 is the pathloss exponent, which 
signifies how fast the channel strength decays as a function of 
distance. 

The pathloss function above is somewhat artificial in its 
dependence on the distances between the receiver user and 
the transmitting base. Technically, it is strictly correct only 
when each user is located vertically to the line connecting all 
bases (as in any horizontal line of Fig.[TJ, which is one possible 
geometry for one dimensional systems |8|. Nevertheless, it has 
the right behavior for — nij as well as for |mi— nij |£ ^ 5. 
The reason we chose to use this model will become apparent 
later. 
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Fig. 1. Schematic figure of wireless network The circles correspond to 
base-stations, while the squares receiver users. The opaque squares represent 
base-user pair that have been "erased" (inactive). 



The users are assumed to try to connect with their own 
nearby cells. The SINR for each connection is given by 



SINRk 



Pkgkk 



'Ej^k 9jkPj 



(3) 



where pk is the transmitted power from the k base to the 
user k and n is the thermal noise, assumed equal for all for 

'in fact, we believe this accuracy of this approximation is better than that 
of a pathloss exponent a > 2 as a free variable, whose value comes out of 
curve fitting of large amounts of data with sizable eiTors | 9j. 



simplicity. For the connection to be possible, a minimum value 
of the SINR has to be attainable i.e. 



SINRk > 70 



(4) 



Therefore, each base should adjust its power to meet this 
criterion. As a result, the following set of equations should 
be simultaneously met 



n ^7o ^9kkPk 



n ^gjkPj > 1 



(5) 



The above equations constitute a set of linear (planar) con- 
straints on the powers. The minimum total power is reached 
at the apex of the conical section of allowed powers, given by 
equality of all constraints above. Defining the matrix M with 



elements Af^ as 



M,;, 



i=j 



we may write the above equation as 



(6) 



(7) 



where P is the vector of powers that satisfies the equality 
constraint above and J = [1,1,...,!]-^. As a result, the 
minimum average power is given by 



Pa 



N 



(8) 



Clearly, for the inverse M to be well-defined, all its eigenval- 
ues have to be positive. 

Now, due to the circulant structure of the matrix, its eigen- 
vectors will be equal to Vq, which at position m has the value 



(9) 



where the vector q resides in the fundamental cell of the 
reciprocal lattice of the lattice m 



2tt 



(10) 



with < kg < L flOl. The eigenvalue corresponding to the 
vector q is simply the Fourier transform of any line of the 
matrix M 



A(q) 



7o 



m- 



(11) 



It is important to note that these eigenvalues are real because 
the matrix M is real and symmetric, which is the reason for 
choosing ([T]i. Finally, since the vector J is proportional to the 
q = eigenvector of M we conclude that 

1 



Pa 



"70 



Ao 1 - 70 I] 



(12) 



m^^O (|m|2+s2) = /2 



where Ao = A(q = 0). This solution corresponds to all bases 
transmitting with the same power. The condition for finite total 
power in the case of no erasures is simply 



70 



^ (|m|2^2 + ^2)a/S 



< 1 



(13) 



III. Sparse Network 



The above equations are valid in the absence of erasures 
of base-stations. Erasures may take place in two cases. First, 
perhaps in a particular location a base station may not have 
been erected. Therefore, no mobile is expected to receive any 
data in that cell. Second, in a cell with a few users present, 
once in a while all data buffers for all connected users may be 
empty. Hence the corresponding base-station will be silent. In 
principle, it may still be transmitting a pilot signal, in which 
case the required SINR will be much smaller. For simplicity, 
we will assume that the connection to its cell-user(s) will be 
off. As a result, we are faced with a situation, in which, each 
row and column, with probability e has zero entries. Thus the 
matrix we now need to analyze is 



Mr ^ EME 



(14) 



where E is a diagonal random matrix with independent 
diagonal entries Ci on the ith entry, with 



P[e, = 0] = e 
Ph - 1] = 1 - e 



(15) 



The corresponding equation for the power vector becomes 

EMEP = EJ (16) 

where we have made sure to also erase the appropriate entries 
from J on the left-hand-side by multiplying with E. However, 
simply inverting the matrix Mc is problematic, since it has 
approximately A^c zero columns and rows and corresponding 
zero eigenvalues. Instead, we regularize the matrix by adding 
a small positive constant on the diagonal, namely 



P = [EME + eI]"^EJ 



(17) 



To obtain the average power as before we multiply on the left 
with J^^E to collect all the components with non-zero power. 
Thus we have 



1 



1 



Pav 



lim — J^E [EME + el] ^ EJ 



(18) 



We will state the result below, and defer the proof to Appendix 

El 

Theorem 1 (Average Minimum Power). Let (3 be the solution 
of the equation 



(271)'^ /3 + A(q) 



(19) 



where A(q) is the q eigenvalue of the matrix bM ( 1771 ) and the 
integral is over the fundamental cell of the reciprocal lattice of 



the system. Then the average power per base-station is given 
by the equation 

1 1 



Pa 



(20) 



1 - e /3 + Ao 

Remark 1.1. In the appendix the above theorem is proved 
straightforwardly for one-dimensional lattices using results 
from |6|. We will argue that this is valid for two-dimensional 
lattices and for all feasible 70 through simulations and defer 
the complete proof for a longer version of the work. 

It is easy to see that the above result reduces to (fT2l l in the 
limit that e — ?> 0. Indeed, when Ao > 0, the right-hand-side of 
( fT9l l is zero only when /3 = 0. Hence (fT2l) follows. We can 
also show that for e 1 the average power simply takes its 
non-interference value 
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(21) 



In the presence of erasures, the maximum SINR target that 
is supported in the system is reduced to 



7o 



70 



(22) 



1 + n/3c7o 

where /3c is the solution of ( fT9l ) and 70 is given by the equality 
of O. 

In addition to the mean power, we can also obtain the 
variance of the powers in the network. Again the proof is 
deferred for Appendix [A] 

Theorem 2 (Variance of Minimum Power). The variance of 
the powers of the individual bases in the network is given by 

f '^q 

VariP)^-^^-^LS2^P^^^ (23) 
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A(q) 



(2^)" (/3+A(q))^ 

A. Analysis of one dimensional array 

We will start by analyzing the results for the simpler one- 
dimensional case, where the base-stations are located in a ring. 
The numerical results are plotted in Figures |2l [3] 

First, we see that the average power obviously is a increas- 
ing function of 70 and a decreasing function of c. 

Second, we see, both analytically and numerically, that the 
average power is finite for a range of 70 above the value for 
which Ao ~ 0, which is the value of 70 for which the c = 
solution diverges. However, this comes at a cost: For such 70 
for which Ao < a second solution to the fixed point equation 
appears (plotted, for the case of a = 4 in Fig.|3]in green). This 
value of Pave is unstable however, having a negative variance 
( |23] l. This second solution of pave persists until it merges with 
the stable value at a critical value of 70. However, at that point 
the variance of the stable solution diverges. Beyond this point 
no finite average power can be supported for the given value 
of e. Thus we see a discontinuous transition at that point, 
from a finite average minimum power to an infinite value. 
This is because beyond this point, at least one power has to 
be infinite. In fact, the finite value of pave all the way to 
the critical value suggests that at most a vanishing fraction of 
bases have unbounded powers. 



Average power vs target SINR; PL =4; noise=1 ; d=0.5 




Fig. 2. Plot of average power as a function of the target SINR for various 
values of e. The red curves correspond to the numerically generated values. 
The path-loss exponent is set to « = 4, while the noise level n = 1 and the 
ratio s = 5/1 = 0.5 



PL =4; n=1:d=0.5 




Fig. 3. Plot of average power as a function of the target SINR for various 
values of e for the parameter values of Fig. [2] In addition to the curves plotted 
also in Fig.[2](blue) we plot the curves corresponding to the unstable solutions 
of )19t (green) and the variance of the power for the stable solution (red). 



Another interesting result from the numerical analysis is that 
beyond the value of 7 where Ao = the behavior becomes 
sample dependent. Thus the value at which the power becomes 
infinite varies somewhat from sample to sample (as well as 
the whole curve beyond the value where Ao =0). As we shall 
show in the Appendix, the end point of this curve has to do 
with the largest eigenvalue of the matrix EME. For finite 
but large N it is known for similar random matrix systems 
lim that although the eigenvalues of the matrix converge fast 
to their asymptotic spectrum, the largest eigenvalue fluctuates 
much more outside the support of the asymptotic spectrum. 
Hence, the theoretical curve seems to be the upper limit, 
beyond which no finite power can be supported. The numerical 
curves in Fig. |2] especially for the larger values of e have 
been picked after running several simulations and picking the 




Fig. 4. Plot of average power as a function of the target SINR for various 
values e for a two-dimensional square grid of base-stations. The path-loss 
exponent is a = 5. The curves have the same meaning as in the previous 
figure. 



one with the largest critical value. How can we explain this 
behavior? For example, if a region happens to have a large 
concentration of bases it may create a situation where the 
power of a particular base diverges. Such a happenstance may 
be less probable if each base is affected by a "mean-field" of 
all bases, including far away ones. 

B. Analysis of two-dimensional grid 

It is remarkable that these results also hold for two- 
dimensional grids, even though we have not provided a proof 
in this short paper version. This can be seen in Fig. |4] where 
the average optimum power for a square grid 50 x 50 of base- 
user pairs is plotted as a function of 7 for different densities of 
erasures e. We plot the analytic curve as well as the numerical 
average over approximately 500 realizations (solid line) and 
the curve for the realization that has lasted to the largest 7 
before becoming unfeasible (dashed). We once again observe 
the same behavior. The "average" curve stops at some point 
after the value of 7 for which Ao(7) = while the dashed 
lines go up to close to the analytic curve. In Fig. |4] we also 
plot for comparison a "mean-field" curve, which assumes that 
the interference is spread evenly over all neighboring sites, 
such that 

P9ii 

7 — — ^ (24) 



It becomes clear that while at small 7 it behaves decently, it 
significantly overshoots the singular point. 

IV. DISCUSSION 

In this paper we have studied the optimal power vector that 
achieves an SINR target criterion in a wireless network setting 
where both randomness and interference are relevant. We have 
applied methods from Random Matrix Theory to obtain the 
average optimum power and its variance. Our numerical results 
show remarkable agreement with predicted analytical calcula- 
tions, both for one- and two-dimensional network geometries. 
We observe that in the presence of randomness the optimal 



power does not go to infinity in a continuous fashion. Rather, 
beyond a certain point, no finite power solution exists. We 
also numerically observe that beyond the point where the pure 
system has finite optimal power, real random systems become 
unstable, making finite power unfeasible for increasing SINR 
targets at different points. The ones that remain finite reach our 
analytical limit. This suggests that randomness creates sample 
to sample fluctuations, which make local clusters of cells and 
hence the whole network unfeasible to operate at that target 
SINR. Therefore, it does not make sense to look for globally 
optimal power vectors beyond a certain network size. Instead 
sacrificing the connection of a few "bad" bases by adding some 
outage may increase the overall behavior dramatically. 

Appendix A 
Proof of Theorems[T]and[2] 

Both proofs will be assume one-dimensional networks. The 
case of d = 2 is deferred to a later publication. 

Proof of Theorem [7} We start by defining the matrix A 

as 

A = EF (25) 

where F is the Fourier matrix that diagonalizes M = Fi^AF 
and denote its column vectors by a(q) and, in particular, 
a(0) — iV~^/^EJ = ao. As a result, ( fTSl ) may be rewritten as 



AoaoaJ + AAqA^ + el 



ao 



(26) 



where Ao is the diagonal matrix of eigenvalues of M with 
the one corresponding to q = (i.e. Aq) set to zero. After the 
application of the matrix inversion lemma we obtain 

1 .. 1 
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(27) 



1 - e el^o+ /?(e) + Aq 
where 

/3(e)-i=aS[AAoAt+el]"'ao (28) 

In Q it is shown that the above quantity takes (a.s.) a 
deterministic value in the large N limit, which is the solution 
of the following equation 

4 H\ ) ^29) 



JIT) ■> ^27r)'' + A(q) 

Taking the limit e ^> 0+ in the above and setting /3(0) = /3 
concludes the proof. ■ 
We now turn to finding the variance of the powers 

Proof of Theorem^ We start by noting that the normal- 
ized second moment of the power vector P is given by 
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lim aj [EME + el] ^ ao 



(30) 



lim — at [EME + ell^^ao 

e^0+ oe 



lim 



1 - c .^0+ (/3(e) + Ao)' 

where /3'(e), the derivative of /3(e) with respect to e can be 
evaluated directly from (|29] l. Subtracting off p^^^ concludes 
the proof. ■ 



It is interesting to point out that the condition that the 
denominator in ( l23T l vanishes, which signifies the end of the 
feasibility of a solution, is equivalent that the largest edge of 
the asymptotic spectrum of the matrix EME is zero and hence 
/3 cannot admit a real solution. To show this we start from the 
realization the fact that M and E are asymptotically free [71 • 
Then defining the Cauchy transform as 



F[z) 



(31) 



z z 

where r\ is the r\ transform of the matrix |[l2l . Taking advantage 
of asymptotic freedom Q, lfT2ll we can show that 

rfq e/f(z) 
(27r)'^ A(q) - z + e/F(z) ^ ^ 

Identifying e/i^(0) as /3 we see that if z = is not part of 
the asymptotic spectrum, /3 has a real solution in the above 
equation, which coincides with ( fT9] l. On the other hand, when 
z = is inside the support of the spectrum, /3 will also have 
an imaginary part. At the edge of the spectrum the condition 
ImF{zedge) = implies 

dq A(q) — z 



(27r)'^(A(q)-z + /3)2 
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